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' Abstract. There exists a large class of quantum many-body systems of Calogero-Suthcr- 

ry | ■ land type where all particles can have different masses and coupling constants and which 

nevertheless are such that one can construct a complete (in a certain sense) set of exact 
eigenfunctions and corresponding eigenvalues, explicitly. Of course there is a catch to this 
result: if one insists on these eigenfunctions to be square integrable then the correspon- 
ding Hamiltonian is necessarily non-hermitean (and thus provides an example of an exactly 
solvable 'PT-symmetric quantum- many body system), and if one insists on the Hamiltonian 
to be hcrmitean then the eigenfunctions are singular and thus not acceptable as quantum 
mechanical eigenfunctions. The standard Calogcro-Suthcrland Hamiltonian is special due to 
. the existence of an integral operator which allows to transform these singular eigenfunctions 

£3 into regular ones. 

O ! 
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1 Prologue 

In summer 2003 Stefan Rauch informed me that a colleague of us would soon visit him in 
^ ■ Linkoping and thus pass through Stockholm, and he suggested that I should invite him for a 

seminar: this colleague was also working on quantum integrable systems of Calogero-Sutherland 
(CS) type [H [2] and I would want to talk to him. I thus got to know Vadim Kuznetsov on 
September 1 in 2003, short before his seminar where he explained to us their work on the Q- 
operator and separation of variables for the CS model [3]. Stefan Rauch was right: not only did 
I find the seminar very inspiring and interesting, but also in the discussion afterward I learned 
about a (for me) new motivation for studying CS type systems which, since then, has influenced 
and helped me significantly in my work on this topic. By 'learned' I mean here that some 
way of looking at things, known only from hear-say before and very different from my own, 
suddenly 'sank in' and made me see things from a new angle. This was a rare experience for me 
which, I felt, was due to a special channel of communication and mutual understanding I had 
when talking to Vadim. Vadim was enthusiastic about our science (not only then but in all the, 
unfortunately too few, instances I had the privilege to meet him), his enthusiasm spread to me, 
and this was a main reason why I soon afterward returned to my work on the CS system [H [5] 
and thus stumbled over results much more extensive and beautiful than what I had expected 
and for which I am very grateful (these results were announced in [61 [7] , and the final version of 
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the detailed paper on this, by which I intend to replace [8], is still in the stage of being 'nearly 
finished' since 20 months. I miss Vadim's interest and kind encouragement, but I intend to 
finish it soon, I want to make it as good as I can, and when eventually finished it will be to a 
large degree Vadim's paper since it would not exist without him). 

In the present paper I will present a story which I had carried with me since quite a while (my 
old draft on it which I use when writing this paper was dated in October 2000) without being 
able to fully make sense of it. It is about a seemingly provocative result, closely related to this 
above-mentioned method to construct explicit formulas for eigenfunctions of such systems [7], 
and I discussed it with Vadim on March 2 (or 3) in 2005, shortly after a seminar to which 
Vadim had invited me to explain them my approach. I brought up this story in response to 
a question of Vadim, and after our discussion I understood its place and it finally made sense 
to me. I remember Vadim liked this story, and this is why I choose it for the present article. 
Moreover, as will become clear later, this story is related in spirit to Backlund transformations 
and the Q-operator which have played a central role in Vadim's research; see e.g. [9} 110] . 

My plan for the rest of this paper is as follows. In the next section I will summarize some 
known facts about CS systems which are relevant to my story. In Section 13.11 1 will present 
the 'strange' exact eigenfunctions for a large class of generalized CS-type systems, and in Sec- 
tions [3?2] and [3T3] I clarify their nature by pointing out a natural interpretation in the context of 
'PT-symmetric mechanics [IT]: they provide non-hermitean "PT-invariant quantum many-body 
Hamiltonians with real spectrum and computable eigenvalues and eigenfunctions. In Section U] 

1 will explain that these strange eigenfunctions for the CS differential operator are interesting 
and useful since they can be transformed into regular eigenfunctions. This provides an alterna- 
tive approach to results previously obtained in [7]. I will end with a few remarks in Section 
Technical details and some more technical proofs are deferred to three appendices. 

2 Background 

2.1 CS model and Jack polynomials 

The Calogero-Sutherland (CS) model [HE] is defined by the differential operator 



with Xj € [— 7T, 7r] coordinates on a circle, N = 2,3, . . ., and A > 0. This differential operator 
essentially defines a quantum mechanical model of N identical particles moving on a circle and 
interacting with the 1/sin 2 two-body potential, and the parameter A determines the coupling 
strength. The CS model is famous among theoretical physicists since it can be solved exactly: 
the differential operator Hjy has exact eigenfunctions of the form 



which are labeled by partitions n of length N, i.e., n = («i,n2, ■ ■ ■ , «jv) with rij integers such 
that 

n\ > ri^ > • • • > nj\r > 0, 

N 

and where P n are certain symmetric polynomials of degree |n| = ^ rij in the variables 




(1) 



tf n (x) = tf (x)P„(z) 



(2) 




(3) 
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with 




the groundstate eigenfunction. Moreover, the corresponding exact eigenvalues have the following 
remarkably simple form, 

N ( 1 \ 2 
E n = Yl [n j + ^N + l-2j)\\ . (5) 

The functions -P n ( z ) are called Jack polynomials and can be characterized, somewhat informally, 
as follows 




+ lower order terms 



where the sum is over all distinct permutations P of (ni,7i2, • • • ,nj\i); see e.g. |12] for a precise 
characterization. 

The following two remarks are on technicalities which can be often safely ignored, but I regard 
them useful as preparation for our discussion of the generalized CS models in the next section. 

Remark 1. For non- integer A our definition of \&o(x) above is only complete in the wedge 
—it < x\ < X2 < • • • < xn < tt where we require it to be real (since the phase factors (— 1) A 
obtained by permuting the arguments are obviously ambiguous). There is one natural method 
to fix this ambiguity by analytical continuation: extend the definition of l I / o( x ) to the other 
regions by continuing the Xj to the complex plane, Xj — > Xj + iej with e > 0, performing the 
necessary permutation, and then taking the limit e — > 0. 

Remark 2. It is worth mentioning another seemingly technical point which, however, will 
play an important role for us later: our 'definition' of the CS model above by Hn in (pQ) was 
somewhat vague since, to be precise, this quantum mechanical model is defined by a particular 
self-adjoint operator on the Hilbert space L 2 ([—tt,7t] n ). Such an operator is only uniquely 
defined if one specifies an operator domain, and the latter we did implicitly by characterizing 
the eigenfunctions. To see that this is not only a technicality we note that the CS model is 
specified by the parameter A, but the coupling constant 

7 = 2A(A - 1) 

in (HJ) is invariant under A — > 1 — A. Thus in the coupling regime — 1/2<7<0 one and the 
same CS differential operator corresponds to two different CS models. In fact, the differential 
operator in ([1]) has many more self-adjoint extensions. We will later encounter other examples 
where one and the same differential operator defines different Hilbert space operators. 

We note that these Jack polynomials are just one example in a whole zoo of symmetrical 
polynomials related to eigenfunctions of CS-type models and which naturally generalize the 
classical orthogonal polynomials (like Hermite, Laguerre, Jacobi, . . . ) to the many-variable 
case. These polynomials were, and still are, extensively studied in mathematics, and many 
beautiful results have been discovered; see e.g. [13] and references therein. It was only due to 
discussions with Vadim that I became aware that my results could be relevant in that context, 
and it was he who convinced me to write [7] where I could obtain explicit formulas for the Jack 
polynomials which (to my knowledge) were not know before. One aim of the present paper is to 
give an alternative derivation and interpretation of these formulas. To be specific I restrict my 
discussion to the Jack polynomials and only mention in passing that all results allow for a (rather) 
straightforward generalization to the other symmetric polynomials mentioned above [141 115] . 
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2.2 A generalized CS model 

There exists a generalization of the CS model describing distinguishable particles and which 
is partially solvable in the sense that is has a groundstate eigenfunctions and corresponding 
groundstate energy which can be computed explicitly. This model is defined by the differential 
operator 

j =1 J i<fc 4sm I— (Xj — Xk) 

with arbitrary mass parameters Mj/2 > and the coupling constants 

ljk = (Mj + M h )X(MjM k X - 1), (7) 
and it has the following exact groundstate 

$o(x) = l[ S ml-(x k -x j )) (8) 

j<k v / 

with corresponding groundstate energy 

3 N 

(9) 





Hiv^o(x) =5 *o(x). (10) 

This fact is known since quite a while (see e.g. [E]), but due to its importance for us a short 
proof is given in Appendix [A] This fact plays a twofold role in our story: firstly, it suggests 
that the Hamiltonian in §§§ might have further exactly computable eigenstates, and the search 
for these led to the result presented in the next section, and secondly, it implies, as special case, 
the following result which will play an important role for us in Section S) 

Lemma 1. The function 

" Tl sm((l/2)( Xj -x k )) x n sin((l/2)( yi -y fe )) A 

rn / \ ^ 2-/ W yj)\<j<k<N l<j<k<N 
F N (x;y)=ce 



JV 

n ^((1/2)^-^))^ 
j,k=i 

with Xj, t/j complex variables, c and P arbitrary constants, and the CS Hamiltonian = i2jv(x) 
in ([1]), obey the following identity, 

H N ( X )F N (x,y) = H N (y)F N (x,y). (11) 

(The proof is given in Appendix lB.il ) 



As was pointed out to me by Vadim, the identity in (jlip is also implied by a well-known 
generating function of the Jack polynomials which is closely related to the functions F^', see 
Proposition 2.1 in [12]. We will use this identity to construct an integral transform Fn in 
Section U] which commutes with the CS Hamiltonian and which, for this very reason, will be 
useful for us. This integral transform is similar to the Q-operator defined and exploited in [3], 
and the usefulness of defining this integral transform Fn was suggested to me by Vadim. 
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3 Strange exact eigenf unctions of CS-type system 

In this section we construct and discuss singular eigenfunctions of the CS-type differential opera- 
tor in ([6|) . We note that this result highlights the importance of square integrability in quantum 
mechanics: if one ignores this condition it is actually easy to construct explicit eigenfunctions for 
a large variety of models. We first will discuss the simplest class of eigenfunction (Section l3.ip . 
Section 13.21 contains an interpretation of these as eigenfunctions of a non-hermitean Hilbert 
space operator with purely real spectrum. We also sketch a generalization of this construction 
illustrating that one and the same differential operator can define a large class of different 
Hilbert space operators. 

3.1 Construction of singular eigenfunction 

We now construct exact eigenfunctions of the differential operator in ([6]) with arbitrary coupling 
constants jjk- These eigenfunctions are labeled by integer vectors n £ "L N and are linear 
combinations of the monomials 

/ n (z) := zfzf ...z$, n+ = n . + Sj (12) 

where s £ is arbitrary for now. It is important to note that these eigenfunctions are not 
symmetric, i.e., not invariant under permutations of the particles. Moreover, we will find that 
these eigenfunctions are only well-defined in a region contained in the following domain, 

Q N = {z e C N \ \zi\ < 1 22 1 < • • • < |zat|} 

and thus are, in particular, not elements in the Hilbert space L 2 ([— 7r, 7t]^)0 To state our result 
we will also need the following partial ordering of integer vectors 

m ■< n 44> rrij + m^+i + • • • + mjy < fij + n-j+i + • • • + n^r V j = 1, 2, . . . , N, 

and we will use the special integer vectors 

Ej/j — Gj e^ 

with ej the standard basis vectors in Z , i.e., (ej)e = Sj£ for all j, i = 1,2, . . . , N. We will also 
use the notation 5 n (m) for the Kronecker delta in Z N , i.e., 

<^n(m) = &ni ,mi ^n.2,r«2 ' ' ' ^nff,m,ff 

Moreover, we will use the following subset of Zi N : 



I j<k 




Note that n + fi ^ n for all p, G C_ and n £ 

We start our construction by the following simple observation: 

Lemma 2. In the region H,n the differential operator in ([6]), for arbitrary jj^, acts on the 
functions in (|12p as follows, 

oo 

HnU = E n f n - ^ Ijk ^2 v fn+"V ]k (13) 

j<k v=l 

where 



1 In Section [3^2l we will discuss an alternative interpretation of these functions as elements in this Hilbert space. 
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Proof. Using 



, d _ d ^ 1 Zj/z k __V^ Z/ 

dxj 3 dzj 4sin 2 (l/2(x i - x k )) (l-Zj/z k ) 2 ^ \z k/ 

for \zj/z k \ < 1, we can write the differential operator in ([6]) on the domain fl^v as follows, 

N 1 / a \ 2 



*» = E^ -EwE"(|) 



J 

This implies the result. 



Equation (113j) shows that the action of Wtv on the functions / n has triangular structure 
in the following sense: Ti^fn is a linear superposition of functions / m with m ^< n. This 
suggests that T~Lm should have eigenfunctions P n which can be expanded as linear combinations 
of functions / m , m ^ n, and with eigenvalues E n . This turns out to be the case, and one finds 
that the corresponding expansion coefficients can be computed by diagonalizing a triangular 
matrix if a certain non-degeneracy condition is fulfilled, and that this matrix is simple enough 
to allow for a fully explicit construction of these eigenfunctions. It is also possible to establish 
a necessary condition for absolute convergence of this series. To summarize: 



Proposition 1. Let E n be as in (|14p . 

N 1 

&n(m) = E m - E n = Jf( m i ~ n i)( m i + n i + 2s i)> ( 15 ) 

j=l 3 

and n G 7L N , s £ 1^ and Mj > such that 

|6 n (n + A)|>A Vfi£C N (16) 
for some A > 0. Then the function 

P n = fn+^2a n (m)f m (17) 

with 



„(m) = 6 n (m) + E II ( E yirkr E v r I - — r-^-r — — r (18) 



s=\r=\ \j r <k r u r =l / [7 & n(n + Yl u ^jik l 

r=l V t=\ 

is an exact eigenfunction of the differential operator in © corresponding to the eigenvalue E n 
in (|14p . Moreover, the infinite series defining this function is absolutely convergent in the region 
where 

EM(T^<- (»> 

Proof. Inserting the ansatz (fTT|) in the eigenvalue equation (Ti^ — E)P n = 0, using (fT3|) . and 
renaming one summation variable we obtain 

E I (E m - E)a n (m) - ^7jfe E ra n(m - i/E jfc ) J/ m = 

m^n \ j<k v=\ / 
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with a n (n) = 1. Since the functions / m are linearly independent we conclude that this is true 
if and only if 

oo 

(E m - E)a n (m) - j jk ^ z/a n (m - uE jk ) = 

j<k v=\ 

for all m -< n. Since a n (m) = for m > n, solving the latter equation amounts to computing 
a particular eigenvector of a triangular matrix indexed by integer vectors, as anticipated (see 
Appendix [Cj. We conclude E = E n and that the coefficients a n (m) can be computed from the 
following recursion relations 



ctn m 



. oo 

77777 E 7j ' fc E ra "( m ~ ^i*)' " n ( n ) = 1 



nV ; j<k w=l 

using the notation in (|15p . provided that 6 n (m) ^ for all m -< n. The latter is guaranteed by 
our assumption in (fT6j) . We can solve this recursion relation by iteration (see Appendix [U]). We 
thus obtain 

oo 



"n( m ) = E"" ( m 



s=0 



where 



ak 0) (m) = <5 n (m), a4 s) (m) = — — - V" j jk V] i/a4 s _1) (m - z^E 

o n (m ^— ' 

nV ; j<k v=\ 

for s = 1, 2, . . .. The formula in (|18p is obtained by computing 



( m ) = jn^s E E 



^s— 1 J j 3 _i</e s _i 



=1 



1 



b n ( m - ^ z^E iA ) j s - 2 <k s -2 
t=s-l 



E ljs-2k B -2 E 



^-2=1 fe n m- X] ^ E iA 



1 



and replacing m in the arguments of the fe n by n + v r^j r k r - 



r=l 



To see where the eigenfunctions P n are well-defined insert (|18p in (|17p and obtain 

1 



s=lr=1 \ jr<kr Ur= i J Y\ fe n ( n + J2 v l^j t k 

r=l V 1=1 

where use the notation / m (z) = z m+ . We can estimate this using (|16p 

l^n(z)| < \Z 



n+ J2 VrEjrkr 
Z V r = l 



(20) 



S = l r = l V j r <k r 



u r =l 



r=l 



Zkr 



5l 









= E EmE^ 

which implies absolute convergence of the series defining P n provided that the condition in (|19|) 
holds true. ■ 
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Remark 3. As the proof above shows, the condition in (|19p is sufficient but far from necessary. 
We believe that, if b n (m) ^ for all m ^ n, then -P n ( z ) is well-defined in all of Qn. 

Remark 4. We call the eigenfunctions in Proposition Q] singular since they seem to diverge as 
Zj — > Zk, at least for positive coupling values: from a heuristic argument we expect that they 
behave like 

„ \ z . _ z ,|(i/2)- v / (T7ij+^ MjM k 

This argument suggests that these eigenfunctions are elements in L 2 {[— 7r,7r] ) if 
all jjk are negative and larger than —1/4. It would be interesting to prove that this was indeed 
the case since then these 'singular' eigenfunctions would be actually acceptable as quantum 
mechanical wave functions. 

Remark 5. The idea that the (standard) CS differential operator has eigenfunctions as in (jTTJ) 
was probably first used by Heckman and Opdam; see e.g. |17j . 

Remark 6. Other remarkable special case of the generalized CS differential operator in © 
and ([7|) allowing for polynomial eigenfunctions are for Mj = 1 for 1 < j < n < N and Mj = — 1/A 
for n < j < N, as discovered in [18] for the case n = N — 1 and explored for general n in [191 120]. 
In fact, all our results on the CS model eigenfunctions described here can be extended to these 
deformed cases 1151. 



The following two subsections can be skipped by readers not interested in "PT-symmetric 
quantum mechanics. 

3.2 Exactly solvable PT-symmetric quantum many-body systems 

We now show that the eigenfunctions in Proposition Q] can be interpreted as solutions of an 
(partially) exactly solvable "PT-symmetric quantum many-body system; see e.g. |llj . For that 
we observe that complex variables (|3|) with \zA = R J , R > 0, can be obtained by shifting the 
particle coordinates to the complex plane as follows, 

xj — > Xj + ije 

and doing this replacement in the differential operator in ([6]) we obtain 

N „ 2 

«»-Efb+E y < 21 > 

j'=i 3 i j<k 4 sin 2 ( -(xj — x k — i[k — j]e) J 

It is interesting to note that this defines an operator on the Hilbert space L 2 ([—tt,it] n ) which 
is not hermitean but rather obeys 

1-/* 'i-i 

n N — n N 

where the star means Hilbert space adjungation and the bar complex conjugation. However, this 
differential operator is invariant under the following combined parity (V) and time reversal (7~) 
transformation [11] : 

""PI . ^ ' ^ ^ ^ 

implying pj = —id/dxj — > pj. Moreover, Proposition [1] provides exact eigenfunctions of this 
differential operators, and the corresponding eigenvalues are indeed real, as they should be [TT] : 
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Corollary 1. The functions in ()20p with 



N 



+ 



E (iasj— J'e)("y+Si) 



Z 



m 



= e J 



=1 



(22) 



and 6 n (m) m ()15p are exaci formaJQ eigenf unctions of the VT -symmetric Hamiltonian TIn 
in ()2ip corresponding to the eigenvalue E n in (|14p. Moreover, this eigenfunction is a well- 
defined function in L 2 ([—7r,7r] N ) provided that (|19|) ZioZds £rwe /or some A > and 



Sketch of proof. The only non-trivial part is to verify the sufficient condition for square in- 
tegrability: compute 



a n (m) in (118p . and majorize by a geometric series, similarly as in the proof of Proposition [T] 



It is interesting to note that for the CS differential operator in (pQ) and Sj in (|27p below, the 
condition in ([23]) holds true automatically for all partitions n and sufficiently large e (e.g. any 
value e > log(-R) with R in (|30p will do), as shown in the next section (see Lemma [3] and its 
proof in Appendix IB. 2h . 

If for such a PT-symmetric Hamiltonian the exact eigenfunction in Corollary Q] exists for all 
integer vectors n, then these eigenfunctions provide a complete basis in the Hilbert space (it 
is not difficult to prove this, using that the functions / n (z), n G Z N , are a complete basis). 
In this case we have an exactly solvable model, otherwise we only might have an partially 
solvable model. However, we believe that the condition in (|19p can be relaxed and that it is 
possible to compute all eigenfunctions explicitly, i.e., the PT-symmetric Hamiltonian in ([6]) with 
eigenfunctions of the form (I17p is always an exactly solvable model. 

Note that we have obtained for the differential operator in (|2ip a family of eigenfunctions 
depending on parameters Sj which we are free to choose. This is similar to what is well- 
known for hermitean operators: a differential operator does not uniquely determine a self-adjoint 
operator, but there is usually a whole family of self-adjoint extensions labeled by continuous 
parameters; see e.g. [21]. It is interesting to note that one can easily construct a much larger 
family of "PT-symmetric Hamiltonians corresponding to the differential operator in (|2ip . In the 
following section we outline this construction, but our discussion is somewhat sketchy. We only 
include it here in the hope that it might be of interest in the context of PT-symmetric quantum 
mechanics [11] , 

3.3 Generalized PT-symmetric quantum many-body system 

This section contains a sketch how to construct eigenfunctions the the generalized CS differential 
operator in ([6]) depending on N(N + l)/2 parameters Sj and Aj&, 1 < j < k < N (the eigen- 
functions in Proposition ([I]) correspond to the special case Xjk = 0). 




2{k-j)e 



(23) 




above. 



Define 




j<k 



(24) 



2 The qualifier 'formal' here means that questions of convergence are ignored. 
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with arbitrary real parameters Xjk- Then the ansatz 
*(x) = *o(x)P(z) 

in the eigenvalue equation TL^ = E^> is equivalent to 7~L' N P = E'P with 
N 1 8 2 ( V, 



i M ^ dx ) jtk V4sm 2 ((l/2)(x, -x k - i[k - j]e)) 



Mj + M fc 
7 ife = 7ife- M . Mfe A ifc (A; 



A' 



and E' = E — Eq with 



j<k<£ 3 j^k 3 

Note that for the special case where Xjk = XMjMk, ^jk in (|7|), e = 0, the function ^o( x ) i n (EH) 
becomes identical to the one in flSJ), 7^ fc = 0, and Eq is identical to £$ in ([9]). We thus we 
recover (fTUj) . Another interesting case is Mj = 1, Xjk = X, jjk = 2X(X — 1) in which case Ti' N 
becomes equal to the reduced Hamiltonian used by Sutherland in his solution of the CS model [2]. 
By expanding the cot- and 1/ sin 2 -terms in power series one can compute the action of TL' N on 
the functions in (|12p which, again, is triangular, and one then can construct eigenfunctions 
of the form (|17[) . similar as in Section 13.11 One thus finds the following generalization of 
Proposition ([T|): 

Proposition 2. The differential operator in ^ has formal eigenfunctions labeled by n £ 1, N 
and corresponding to the eigenvalues 



N 



x - {nj + SjY x - ( (rij + Sj ) (w fc + gfc) \ 
E «-2^ M +l^ X M-^f Mk~' + ° 

3=1 3 j<k 3 ft 



with E in (|2oj) . These formal eigenfunctions are given by 

* n (x) = * (x)P n (z) 

iwii/i x I / o( x ) (p4|) and P n ( z ) mi (fT7|) where the coefficients are determined by the following 
recursion relations, 

(S, - S. K(m) = E E fa. ( ^ m' + ^ - ( ""° M " ^ ) - < 26 > 
j<& ^=1 V V i k J J 

with a n (n) = 1. 

Note that the recursion relations in ()26[) still have triangular form and thus can, in principle, 
be solved by the method explained in Appendix O provided there are no degeneracies, i.e. 
E n E m for all pertinent m ^ n. However, the resulting formula will be much more involved 
and probably not very illuminating. 
It is interesting to note that, if 
ry+^j nk + sk 
Mj M k 

then the eigenvalues in ([5]) can be written in the following simple form, 

N P 2 11 

E " = Yl ]jp Pi = 713 + 2 ^ X3k ~ 2 2 v 

3=1 J fe<i k>j 



On Calogero-Sutherland Type Systems 



11 



4 A method to solve the CS model 

In our discussion of the singular eigenfunctions all generalized CS differential operators in §6§ 
could be treated on equal footing. It thus seems natural to ask: is there anything special about 
the CS differential operator in ([1]) with regard to these singular eigenfunctions? To answer this 
question we compare the eigenvalues in ([5]) and (fll]) and observe that, if we choose 

8j = ±(N + l-2j)\ (27) 

and restrict the integer vectors n to partitions, then the singular eigenfunctions of Hn in (|17p 
not only are labeled by the same quantum numbers but also have the very same eigenvalues as 
its regular eigenfunctions in (|2j). Moreover, these singular eigenfunctions are well-defined, i.e., 
the condition in (I16|) is automatically fulfilled. We now show that this is no coincidence: The CS 
differential operator Hjy is special since there exists an integral operator which commutes with 
it and which transforms its singular eigenfunctions in (|17p as specified above into regular ones as 
in §2$. Thus we can obtain all the regular eigenfunctions of the CS model by first constructing 
singular eigenfunctions (which can be done explicitly), and then transforming them using the 
operator Fn. As shown below, in this way we recover the explicit formulas for Jack polynomials 
derived previously in [7J. 

In the following we make these statements more precise. For that we first define the class of 
functions TZjy and Sn to which the regular- and singular eigenfunctions of the CS differential 
operator belongs. 

Definition 1. The regular domain 1Z n of the CS differential operator in ([1]) is the vector space 
of all functions of the form ^q(x)P(z) with x I / o( x ) hi @ and P(z) a symmetric function in the 
variables Zj = e 1Xj which is analytic and bounded on the domain \zj\ = 1 Vj. 

Definition 2. The singular domain Sn of the CS differential operator in ([1]) is the vector space 
of all Laurent series of the form 

E a A/n +A (z) (28) 



for some fixed n £ Z N , with / m in (|12p . Sj in (|27p . and coefficients ap, such that 

N 

i A | Yl R-^ k{k ~ j) < oo (29) 



for some R > 1 which is sufficiently large. 

Remark 7. The convergence condition in (|29p is chosen such that the series in (I28p is absolutely 
convergent in the region where \zj\ = R J Vj, and it then obviously is convergent for \zj\ > R 3 vj. 
We introduce this parameter R > 1 to avoid technicalities: our crude arguments in Appendix lB.2l 
show that the following value of R is sufficiently large, 

R = max(2,N(N -l)\X-l\/8), (30) 

but we suspect that any value R > 1 would do. 

We now summarize the properties of the singular eigenfunctions of the CS differential operator 
following from Proposition [TJ As mentioned, the non-trivial part is that the conditions in (|19p 
are automatically fulfilled. 
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Lemma 3. For all partitions n of length N and Sj as in (|2T[) the series P n in ()17p and (|18p 
/ or Ijk = 2A(A — 1) and Mj = 1 (independent of j, k) are well-defined singular eigenf unctions of 
the CS differential operator ([T|) corresponding to the eigenvalue 

Pn G 5at and H^Pa = E n P n . 

(The proof is given in Appendix IB, 21 ) 

We now explain how to construct the operator F^. The idea is to use the function i*jv(x, y) 
in (jll|) as integral kernel and define 



C 



F N (P)(x) = j ^F(x,y)P(0, % = & (31) 



for suitable choices of the constants P and c, and the integration domain C. Then Lemma [J 
guarantees that this operator commutes with , and it maps eigenfunctions of Hn into eigen- 
functions of Hn without changing the eigenvalues, even though the character of these eigenfunc- 
tions is very different. A straightforward computation (deferred to Appendix IB. 3|) then shows 
how c, P, and C need to be chosen to get a simple and meaningful result, and by inspection one 
finds the domain and range of this operator. To be precise: 



Lemma 4. The prescription 



n u-£i/6c) A N 



u;zf/i Sj in (|27p and i? > 1 sufficiently large, defines an operator commuting with the CS differen- 
tial operator (pQ) and mapping its singular- to its regular domain: 

F N (H N P)(x) = H N F N (P)(x) 

and 

F N (P) G TZ N for all P £S N . 
(The proof is given in Appendix IB. 31 ) 

Combining Lemmas [3] and E] above we can conclude that the functions i ? Ar(P n )(x) are regular 
eigenfunctions of the CS differential operator to the eigenvalues in ([5]), and they are of the 
form 

F N (P n )(x) = * (x)P„(z). 

N 

One can prove that -P n ( z ) is actually a symmetric polynomial of degree n j an d that -P n ( z ) 

is equal to the corresponding Jack polynomial, up to normalization; see [5] and [15] for details. 
Interchanging summation and the integral transform we thus obtain the following representation 
of the Jack polynomials as a linear superposition of the special symmetric polynomials f n = 
F N (fn): 
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Proposition 3. The Jack polynomial labeled by the partition n is proportional to 
P n (z) = «n(m)/ n (z) 



where 



/ m (z) = n [{ 

i=l V* 



n / , M x n (i - 6/6 

"6 (■"' J l<j<k<N 



--Ri 27ri 6 



A 7 " 

n (i - 



r=l 



s 



a n (m)=5 n (m) + ^[2A(A-l)] s n £ E 

s=l r=l \><fc r ty=l / f] 6 n ( n + ^ viEj t k 

r =i v ^=i 

and 6 n (m) = E n — E m with E n in ([5]). 

Remark 8. A non-trivial point in the result above is possible degeneracies: if the partition 
n is such that there is no other partition m such that m j = n j anc l -^n = E m then 

3 J 

the corresponding eigenfunction of the CS model is (essentially) unique and the symmetric 
polynomials thus only can be equal to the corresponding Jack polynomial. However, if this is 
not the case an additional argument is needed. 

Remark 9. It is interesting to note that, in the special case A = 1, the result above gives the 
following integral representation of the Schur polynomials, 



At \ n (1-6/6) 

"6 t U 3 I X <3< k < N 



27Ti^^ / N 

\=& 7 n (i-^/6) 

j,k=l 

this is a simple consequence of results in [22] (see [E]); see also Appendix B in [23]. 

Remark 10. Integral representations of the Jack polynomials involving integral operators some- 
what similar to ours were previously obtained in [24\ [25] . 



5 Epilogue 

Yuri Suris at a meeting in Rome in May 2001 pointed out to me that the (elliptic generalization) 
of the identity in Lemma Q] looks very much like the Backlund transformation for the classical 
Calogero-Sutherland system, but I did not absorb this comment at that time. In January 2006 in 
Rome again it eventually 'sank in' and I eventually learned about Stefan Rauch's pioneering work 
on this [26]. I thus understand now that it was not a coincidence that it was through Stefan 
that I got to know Vadim who has played a central role in developing the idea of quantum 
Backlund transformations (see e.g. [9]) which, in the special case of the CS model, is based 
on the very same identity. I stumbled over this identity by coincidence when working on a 
project where I studied a quantum field theory model (initially with Alan Carey) with the aim 
to learn more about the so-called fractional quantum Hall effect [4], and through the positive 
feedback and encouragement I received from the integrable system's community I started to 
explore its consequences, and now I find myself working more and more on problems related to 
the mathematical theory of special functions. This happened to a large degree due to Vadim. 
He has encouraged me and given me a few glimpses on some of his visions for which I am very 
grateful and which will probably continue to have a strong influence on my work. There will be 
probably more papers in the future where Vadim's name should be mentioned . . . 
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A Groundstate of the generalized CS model 

As discussed in the text, the Schrodinger operator in © is a partially solvable system in the 
sense that its exact groundstate and corresponding groundstate energy can be found explicitly. 
In this appendix we give a precise formulation and proof of this fact. 

Lemma 5. Let <&o( x ) be as in (jSJ) and 

>±_ i 9 ■ „ „ _ dlog$ (x) 



Of = ±7^— + Vj with Vj 



"i dxj 3 3 dxj 

for arbitrary complex parameters A and Mj . Then 

N 1 

Ew q ^ = Hn - £o 

3=1 3 

with Hn the differential operator in (|SJ) and So the constant in §§§ , and this implies (|1U|) . In 
particular, if all Mj are positive then TLn defines a self-adjoint operator on the Hilbert space 
L 2 ([— 7T, tt] n ) bounded from below by So (Friedrichs extension; see e.g. |21] ), and<&o is its ground- 
state. 



Proof. We compute 

Vj = ^2MjM k \(j)(xj - x k ) with 00) := - cot 



and thus, using 



0'(z) = ---d ) (z) 2 = --sm- 2 ' 



4 rw 4 V 2 , 
we obtain by straightforward computations 

with "Htv i n © and the reminder terms 
(*) = A 2 



EES MjMkM^ixj - x k )4>{ X j - x e ) + \ £ E M i M ^ ) • 

3 k^j£^j,k j k^j / 

The non-trivial part of the result thus is that the three-body terms in (*) add up to a constant. 
To see this we symmetrize the first sum in (*) above as follows, 

-2 ^2 MjM k Mi[4>(xj - x k )(j)(xj - xi) + (p(x k - Xj)4>(x k - x e ) + (p(x e - Xj)4>(x e - x k )]. 
j<k<e 

This allows us to use the following trigonometric identity 

4>(xj - x k )(f>(xj - x e ) + 4>{xk ~ Xj)(f>(x k - Xi) + 4>(x e - Xj)(f>(x e - x k ) = -- 
and obtain 

W = T ( 2 MjM k Mi + ^ MjMkiMj + M k ) j , 

V j<k<£ j<k J 

and by straightforward computations we find (*) = So as in Q. 

Since obviously Qj&o = for all j this proves (fTUj) . Moreover, Q^" is the hermitean conjugate 
of and if all Mj > then J2(l/Mj)Q~jQj therefore defines obviously a non-negative, self- 

3 

adjoint operator. This implies Wat > So, and thus $o is groundstate of Hat. ■ 



On Calogero-Sutherland Type Systems 



15 



B Proofs 

B.l Proof of Lemma [1] 

The identity in (|10|) holds obviously true for arbitrary real parameters Mj and complex variables 
Xj (see Appendix [A}. 

We double the degrees of freedom and replace A by 2 A. Setting Mj = 1, M^+j = — 1, and 
XN+j = Vj for j = 1, 2, . . . , A we find that H.2N = Hn(x) — -ffjv(y) (note that all interaction 
terms mixing the variables Xj and yj are zero since (Mj + Mjv+jt) = for j, k = 1, 2, . . . , A; we 
also slightly abuse notation here), <J?o(x, y) — -F/v( x >y) for P = 0, and So = 0. Thus the identity 
in CEI implies JTT]) for P = and c = 1. 

We now observe that, due to translation invariance, Hn(x) has the following simple depen- 

N 

dence on the center of mass coordinate X = ^ Xj/N, 

3=1 

d 2 

2M*) = -^ + ... 

where the dots are terms independent of X (they only depend on the relative coordinates 
Xj — X 7\r for A = 1,2,..., A — 1). This implies that the identity in JTT]) is invariant under the 
transformation 

F w (x,y)^ce iP '( X - y )F w (x,y) 

for arbitrary constants P' = PN and c (the invariance under multiplication with c is trivial, of 
course). 

B.2 Proof of Lemma H 

We observe that 

N / N 



M n + E Hjk'Ejk) = ^2 Hjk[nj - n k + (k - j)X] + 

i< k j=i \ k=j+l 



k<j k>j 

which shows that (I16|) holds true with A = 2A. By a simple computation we find that (|19|) is 
equivalent to 

\ z j/ z k\ „ 1 



^(1-l^-MI) 2 |A-1| 



Replacing \zj\ by i?- 7 we find that (|29p is equivalent to 



and this holds true for sufficiently large values of R > 1, e.g. by the following rough estimates 

r __. Rr 1 R-^-NfN - 1) 
l.h.s. of (33) < — = 1 

j<fe 

assuming i? > 2, we find that the given in ([30]) will do. 
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B.3 Proof of Lemma H 

We express the function F(x., y) in terms of the variables Zj = e 1Xj and £j = e iyj (using sin(a;/2) 
z l / 2 {l — l/z)/(2i)), and by a straightforward computation we find 

iMx,y) = C* (x)(2i)^ JV2 -^( Ar - 1 )/ 2 ) ^^™/2 C -P-(iV + l-2 J )A/ 2 -iVA/2x\ 



\j=l 



n (i-&7&) A N 



x i<i<k<N 



N ll^J 

n a - zj/Ck) x i =i 
j,k=i 

To get rid of an awkward constant and the non-analytical (for non-integer A) factors Zj we 
choose 

c = (2i)W-D/2-^) and P = _^. 



This yields 



n (i-^ fc ) A 



N 

771 / x , T( / n l<j<k<N -r-r A(Af+l-2j)/2 

Fjv(x,y) = *o(x) - [[Cj 

n (i-^/6) A j=i 

From this we see that if we choose the integration paths as 

C: y j = - ije, -it < fa < it for j = 1,2,... ,N 

then £j = RPe^i with R = e e > 0, and thus we find that (|31[) is identical with what is given 
in (|32jh 

It is important to note that i*jv(/ n ) for monomials in (|12|) is well-defined if and only if sj 
is chosen as in (|27p (modulo integers, of course), and Sn is the natural domain for Fn since it 
contains all linear superpositions of these / n (z) which are absolutely convergent on C (recall 
|z m | =Y\Rm on C, and Pty): writing F N (P)(x) = tf„(x)P(z), 

i 

^)=n / ^ mrm 

we can estimate on |^-| = 1, 

n (i - w~ h ) x N 

\P{*)\ < II/.' 

n (i - i =i 

fe=i 



and thus prove that Fn(P)(x.) is in the regular domain IZn provided P is in the singular 
domain Sn. 
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C Explicit diagonalization of triangular matrices 

In this Appendix we explain the method to explicitly diagonalize a non-degenerate triangular 
matrix which we use in the main text. The matrix in the main text is indexed by integer vectors, 
but to simplify our notation here we label the matrix elements by integers J, K € Z, i.e., we 
write the matrix as A = (Ajk) j,kgz- 

We are interested in finding the eigenvector v = (vj)j^z of this matrix A corresponding to 
the eigenvalue E = All. Our assumptions about the matrix are as follows, 

Aj K ^0 only if K h J, Ajj-A ll ^0 VJ >- L. 

The symbol here can mean any partial ordering of the integers. We write ax = Akk for 
the eigenvalues of this matrix. 

The vj are determined by the eigenvalue equation 

^ Ajkvr = Evj 

for J ^ L, where vk = for K y L and vl = 1- For J = L we get E = All = and for 



r Y] A JK v K 



VJ = 

(ar, — a j ) 

which gives well-defined recursion relations due to our non-degeneracy assumption. We can 
solve this by iteration: 

oo 

E(s) 
V J 

s=0 

where 



v 



J - °JLi V .J - , _ x 2^ JK V K V S > 1 



which is well-defined since v^Z ^ = for K >- L. This yields 

oo 1 1 

= #JL + V t r V A JiYl V ^X!iY 2 ■ • • 

U^L-aj)^ (a L -a Kl ) K ^ Ki 

x 7 r z2 A Ks _ lKg 5 K L - 

We thus get the following fully explicit formula for the components of the eigenvector, 

AjkxAr^Ks ■ ■ ■ A Ka _ lKs 5 KsL 



vj 



««+E E 



Note that we keep the sum infinite only for simplicity of notation, but it actually has only 
a finite number of non-zero terms: convergence is no problem here, of course. 
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